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Abstract 



We show that the exact beta-function (3{g) in the continuous 2D model 
possesses the Kramers- Wannier duality symmetry. The duality symmetry 
transformation g = d{g) such that (3{d{g)) = d'{g)(3{g) is constructed and the 
approximate values of g* computed from the duality equation d{g*) = g* are 
shown to agree with the available numerical results. 

The calculation of the beta-function f3{g) for the 2D scalar g^'^ field theory 
based on the strong coupling expansion is developed and the expansion of P{g) 
in powers of g~^ is obtained up to order g~^. 

The numerical values calculated for the renormalized coupling constant 
g'^ are in reasonable good agreement with the best modern estimates recently 
obtained from the high-temperature series expansion and with those known 
from the perturbative four-loop renormalization-group calculations. 

The application of Cardy's theorem for calculating the renormalized 
isothermal coupling constant gc of the 2D Ising model and the related univer- 
sal critical amplitudes is also discussed. 

PACS numbers:05.50.+q, 03.70.+k, 64.60.-i, 75.10.Hk 
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I. INTRODUCTION 



In this paper we study mainly the symmetry properties of the beta- function i3{g) for the 
2D (7$^ theory, regarded as a continuum hmit of the exactly solvable 2D Ising model. In 
contrast to the latter, the 2D g^^ theory is not an integrable quantum field theory. This 
means, in particular, that the theory does not possess the factorized scattering matrix, and 
therefore that the thermodynamic Bethe ansatz method cannot be applied at all. 

Thus, despite the fact that the 2D Ising model at h = can be solved by many different 
methods (see ^ for an excellent review), the beta-function j3{g) of its continuum limit is 
to date known only in the four-loop approximation within the framework of conventional 
perturbation theory at fixed dimension d = 2 Calculations of beta-functions are of 

great interest in statistical mechanics and quantum field theory. The beta-function contains 
the essential information on the renormalized coupling constant g*, this being important for 
constructing the equation of state of the 2D Ising model - for example - which remains still 
a challenging problem, rich in applications. This and other considerations do not allow us 
to regard the 2D Ising model as having fully been solved. 

The 2D Ising model and some other lattice spin models are known to possess the remark- 
able Kramers- Wannier(KW) duality symmetry, playing an important role both in statistical 
mechanics and in quantum field theory . The self-duality of the isotropic 2D Ising model 
means that there exists an exact mapping between the high-T and low-T expansions of the 
partition function In the transfer- matrix language this implies that the transfer-matrix 
of the model under discussion is covariant under the duality transformation. If we assume 
that the critical point is unique, the KW self-duality would yield the exact Curie tempera- 
ture of the model. This holds for a large set of lattice spin models including systems with 
quenched disorder (for a review see 

Over twenty years ago the KW self-duality was shown to be equivalent to a Fourier 
tranformation in target space 0. Also, it has been recognised long ago that self-duality 
combined with some special algebraic properties of a model leads to the existence of an 



infinite set of conserved charges Duality is thus known to impose some important 

constraints on the exact beta- function |TT]|T^. 

The other main purpose of this paper is to develop a strong coupling expansion for the 
calculation of the beta-function of the 2D scalar g^'^ theory as an alternative approach to 
standard perturbation theory. It will then be of interest to match this expansion with the 
results of a four-loop approximation (where possible) by constructing a smooth interpolation 
with respect to g. It is in fact well known from quantum field theory and statistical mechanics 
that any strong coupling expansion is closely connected with a suitable high-temperature 
(HT) series expansion for a lattice model From the field-theoretical point of view 

the HT series are nothing but strong coupling expansions for field models, the lattice being 
considered as a technical device to define cutoff-regularised field theories. 

Recently, the high-temperature (HT) series expansions and perturbative calculations for 
the (7$^ field theory at fixed dimensions d < 4 have been a topic of intense studies (for 
references see below). Computing critical exponents and various critical amplitude ratios 
from series expansion data has a long history going back to the early 1960s. Nowadays there 
are a good number of papers containing a large body of information for the A^-vector model 



defined on different lattices for ci = 2,3,4, 5 and arbitrary |13-15|. It is remarkable that 



2 



the HT series data for the zero field susceptibihty x ^i^nd the second correlation moment /i2 
of the A^-component classical Heisenberg ferromagnet have been extended up to the order 
K'^^{K = J/T), the data for the second field derivative of the susceptibility (X4) being 
available through to the order K^^ . Having been equiped with this information, one may 
try to employ different techniques of resummation of the existing HT series expansions, 
like Pade approximants or more subtle approaches, for computing critical exponents and 
universal critical amplitude ratios pB|JTB| -pT|. It is worth noting that the strong-coupling 
behavior of the (7$^ theory has recently been treated within the framework of a variational 
perturbative approach p2| . 

The paper is organized as follows. In Sect. 11 we set up basic notations and define the 
duality symmetry transformation g = d{g). Then it is proved that (3{d{g)) = d'{g)f3{g). An 
approximate expression for d{g) providing good estimates for g^ (the renormalised fixed- 
point coupling constant along the isochore line) is found. In Sect. Ill the HT series expansion 
data are used to obtain the strong coupling expansion of f3{g) for the 2D 0(iV)-symmetric g^^ 
theory in powers of 1/(7 up to the order g^^. Some numerical estimates for the renormalized 
coupling constant (7^ above Tc are obtained. We then compare the fixed point values found 
to those already known from the four-loop renormalization-group (RG) calculations and 
from the HT series expansions. In Sect. IV we also discuss the application of Cardy's 
formula both for the exact calculation of the renormalized isothermal coupling constant g* 
at Tc and, for some universal critical amplitudes, along the isothermal critical line. Sect. V 
finally contains some concluding remarks. The Appendix presents a simple derivation of the 
correlation length ^ and of the exact beta-function f3ising{T) for the lattice 2D Ising model, 
where the temperature T plays the role of an effective coupling constant, and we discuss 
some of their properties. 



II. DUALITY SYMMETRY OF THE BETA-FUNCTION 

We begin by considering the classical Hamiltonian of the 2D Ising model (in the absence 
of an external magnetic field), defined on a square lattice with periodic boundary conditions; 
as usual: 

H = -J J2 ^i^j (2-1) 

<i,j> 

where < i,j > indicates that the summation is over all nearest-neighboring sites; o"j = ±1 are 
spin variables and J is a spin coupling. The standard definition of the spin-pair correlation 
function reads: 

G{R) =< anao > (2.2) 

where < ... > stands for a thermal average. 

The correlation length may be defined in many different ways, all definitions being equiv- 
alent to each other in the close vicinity of the critical point |T3]. This, in fact, refiects the 



arbitrariness somewhat inherent in any renormalization scheme. The statistical mechanics 
definition of the correlation length is given by [^] 



3 



2 d\nG{p) 

e -^^lp=o (2-3) 

The quantity is known to be conveniently expressed in terms of the spherical moments of 
the spin correlation function itself, namely 

/i, = 5:(i?/a)'G(R) (2.4) 

R 

with a being some lattice spacing. It is easy to see that 

= ^ (2.5) 

where d is the spatial dimension (in our case d = 2). It should be mentioned that the above 
definition of ^ differs from the one used in other related approaches, e.g. |Tl] . 



The 2D Ising model near Tc is known to be equivalent to the g^'^ theory with a one- 
component real order parameter. In order to extend the KW duality symmetry to the 
continuous field theory we have need for a " lattice" model definition of the coupling constant 
g, equivalent to the conventional one exploited in the RG approach. The renormalization 
coupling constant g of the g^"^ theory is closely related to the fourth derivative of the 
"Helmholtz free energy", namely d'^F{T,m)/dm'^, with respect to the order parameter m = 
($). It may be defined as follows (see |1T3| , [14| , ^ and references therein) 



9iT,h) = —^ + 3——^ (2.6) 



where x is the homogeneous magnetic susceptibility 



^ = / d'^xGix) (2.7) 



It is in fact easy to show that g{T, h) in Eq.( p.6| ) is merely the standard four-spin correlation 
function taken at zero external momenta. The renormalized coupling constant of the critical 
theory is defined by the double limit 

g* = \imYim^g{T,h) (2.8) 

and it is well known that these limits do not commute with each other. As a result, g* is a 
path-dependent quantity in the thermodynamic {T,h) plane [1^ . 



Here we are mainly concerned with the coupling constant on the isochore line g{T > 
Tc, /i = 0) in the disordered phase and with its critical value 

gl= g{T,h = 0) = -^4^\,=o (2.9) 

The couplings g'^ are of great interest for calculating the equation of state. Notice that the 
two different fixed points values g^ (defined by the above Eq. (p.9|)) and g'i (defined by the 
analogous limit procedure for T T~) computed above and below the Curie temperature 
differ vastly from each other owing to the broken symmetry of the ordered phase. 
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The "lattice" coupling constant g"^ defined in Eq. (|2.9| ) is in a given correspondence with 
the temperature T^. We shall see that it will be more convenient to deal with a new variable 
s = exp{2K) tanh(i^), where K = J/T. The standard KW duality tranformation is known 
to be as follows 

sinh(2i^) = — (2.10) 
^ ^ sinh(2/s:) ^ ^ 

It follows from the definition that s transforms as 5 = this implies that the correlation 
length of the 2D Ising model (see also the Appendix) = -(jz^ is a self-dual quantity. 
Now, on the one hand, we have the formal relation 

«^^^.(«) (2,n) 

where s{g) is defined as the inverse function of g{s), i.e. g{s{g)) = g and the beta-function 
is given, as usual, by 

i'^-m (2.12) 

On the other hand, it is shown in the Appendix that 

^ds 2s(l — s) , , 



^From Eq.s ( p.ll|) - (|2.13|) , a useful representation of the beta-function in terms of the s{g) 



function thus follows 

^ 2s{g){l - s{g)) 
^^^> (1 + s{g)) {ds{g)/dg) ^ ' ^ 

Let us define the dual coupling constant g and the duality transformation function d{g) as 

^i9) = ^y ~g ^ d{g) = s-\^) (2.15) 

where s~^{x) stands for the inverse function of x = s{g). It is easy to check that a further 
application of the duality map d{g) gives back the original coupling constant, i.e. d{d{g)) = 
g, as it should be. Notice also that the definition of the duality transformation given by Eq. 
( p.l5| ) has a form similiar to the standard KW duality equation, Eq. (|2.10|) . 



It is easy to prove that d'{g*) = ±1. The maps we are looking for have d'{g*) = —1, since 
the opposite sign leads to the trivial solution d{g) = g. This is also shown in the Appendix. 

Consider now the symmetry properties of P{g). We shall see that the KW duality sym- 
metry property, Eq. ( p.lOj ), results in the beta-function being covariant under the operation 
9 d{g): 

(3{d{g)) = d'{g)(3{g) (2.16) 
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To prove it let us evaluate j3{d{g)). Then Eq.( p.l4| ) yields 



ma)) = ''^[IZf.^L (2-17) 

{l + s{g)){ds{g)/dg) 
Bearing in mind Eq. (|2.15|) one is led to 

Pidia)) = , (2-18) 
s{g){\^s{g)){ds{g)ldg) ' 

The derivative in the r.h.s. of Eq. (|2.18|) should be rewritten in terms of s{g) and d{g). It 
may be easily done by applying Eq. ( p.l5|) : 

ds{g) ^ d 1 ^ s'{g) 1 
d~g dgs{g) s\g)d'{g) ^ ' ^ 

Substituting the r.h.s. of Eq. ( |2.19|) into Eq. (|2.18|) one obtains the desired symmetry 
relation, Eq. (|216|) . 

Therefore, the self-duality of the model allows us to determine the fixed point value 
in another way, namely from the duality equation d{g*) = g*. One may now test the 
compatibility of this approach with the standard methods by computing d{g) in the simplest 
approximation. For this purpose let us consider the function s{g) given by Eq. (|A.1CI|) (in 
the Appendix). Making use of a rough approximation, one gets 

/ X 2 24 2 1 2 , , 

s(g) ~ - + — ~ = 2.20 

9 9l-l2/g g-12 ^ ^ 

Combining this Pade-approximant with the definition of d{g), Eq. ( p. 15 ), one is led to 

d{g) = 4^^^ (2.21) 
^-12 ^ ^ 

The fixed point of this function, d{g*) = g*, is easily seen to be g* = 14. As expected, 
d'{g* = 14) = -1. 

Before moving to the next topic, we notice that the above-described approach may be 
regarded as another method for evaluating g*, fully equivalent to the standard beta-function 
method. A systematic way of determining d{g) has not yet been developed, in particular 
because of the specific analytical properties of d{g). Namely, because g = and g = oo are 
not regular points of d{g), this function cannot be expanded in a Taylor series around these 
points. 



III. STRONG COUPLING EXPANSION FOR THE BETA-FUNCTION 

Our next purpose is to develop a strong coupling expansion for the 2D scalar g^"^ field 
theory. It has been mentioned already that nowadays the HT series expansions for g{T > 
Tc,0) are known rather well ||T3|-pO|. Eq.s ( |A.9| ) and (|A.10|) , obtained in the Appendix, 



form the basis of our treatment. These equations relate to s, up to the ninth order in 
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the appropriate expansion parameter. Having been equipped with these formulas, one may 
easily calculate the beta-function (3{g) as a power series in g~^. 

Inserting Eq. ( |A.10| ) into Eq. ( p.l4| ) and performing simple but somewhat cumbersome 



calculations, we are led to the desired asymptotic expansion for I3{g) 

p(g) = -2g + 32 - 64:/g + 512// + 512// - 30720// 

- 172032// + 32768// - 172032/^^ + 32768// + 0(^-^) (3.1) 



^From Eq. (|3.1|) it follows that in the large-^f limit (3{g) —2g + 32, whilst in the weak 
coupling regime one has for g ^ : l3{g) ^ +2g |[^-|^, |25|J2^ . It implies that the continuous 
function P{g) changes sign at least once at some fixed point g*. 

Let us get some numerical estimates for g^ now, from Eq. ( p.lD , and compare these results 
with those found from the HT series expansions and those of the four-loop RG calculations. 
In the standard perturbative approach to quantum field theory at fixed dimension one must 
apply some resummation technique to the expansions of j3{g) and other RG- functions. It is 
interesting that at least in low orders of perturbation theory the l/^f-expansion, Eq. ( p.lD , 
does not require the application of a resummation technique. The most reliable numerical 
estimates of g^ were obtained by means of the straightforward solution of the equation 
/?((?*) = 0, from Eq. ( p.l|) , taken within the (^"^-approximation (without the last two terms 
in g~^ and g~^). The five (and best) subsequent approximation are as follows 

/(I) = 16; g*^'^ = 13.6568; /^^^ = 15.0044; 

g*/^ = 15.0784; /^^^ = 14.7632 (3.2) 

Here the index k in g_^ indicates that k + 1 terms are retained in the fixed point equation 
under discussion. These estimates exhibit a regular behavior, the last value being in very 
good agreement with the most recent estimate g"^ = 14.700 ± 0.017 obtained for the square 
lattice |1I^JT5[| . The estimates obtained after taking into account the g~'^ and g~^ terms 
differ significantly from the above values. This is apparently an indication that 1/g-senes 
also require the application of some resummation technique. 

Another approach to obtain a numerical estimate for g^ is a straightforward solution 
of Eq. ( [A.lOp , given in the Appendix, after setting s = 1. In contrast to the fixed point 
equation, i3{g) = 0, it yields a rather poor value of the renormalized coupling constant, 
g*^ = 12.533, compared to the value reported in P JI^ , p!5[ . 



It is interesting to compare our results with those obtained from the beta-function of the 
2D Ising model and computed in the four-loop approximation, known to provide more or 
less satisfactory results for the critical indices 



(3{v) = 2v- 2v'^ + 1.432346^;^ - 1.861533^;^ + 3.1647764?;^ + 0{v^) (3.3) 



To obtain the beta-function in our normalization we have to change variables |£4 



9 = m = y/5(^) (3.4) 

The analysis based on the Pade-Borel method of resummation of asymptotic series yields 



/ = 15.08 ± 2.5 [Q, which slightly exceeds the best values obtained from the HT series 
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calculations: g*^ = 14.70 ±0.017 gl = 14.67 ±0.04 ^ (it is tempting to conjecture 



24 



that gl = ^ 

The exact numerical value of the critical exponent u = f3'{g*) governing the leading 
corrections to the scaling laws |23,25] is still unknown. The straightforward application of the 
above strong-coupling expansion results, Eq.s ( |3.1| ) and (|3^), yields (within a conventional 
uj>0 definition) u = (3 {g* = 14.76) = 1.88. This estimate was obtained without exploiting 
a resummation procedure. It differs vastly from the estimate found from the four-loop RG 
calculations combined with the Pade-Borel method: uj = 1.3 ± 0.2 p4|. This value agrees 
well with that predicted by conformal field theory: = | 



IV. ISOTHERMAL COUPLING CONSTANT AND CRITICAL AMPLITUDES 

The two preceeding Sections were devoted to computing the approximate value of the 
renormalized coupling constant g^ at h = in the isochore limit. Here we remark that 
in two dimensions there is a possibility of calculating the exact value of the renormalized 
coupling constant g* in the isothermal limit, i.e. at the Curie point in an applied magnetic 
field, namely 

g: = ]img{T = T,,h) (4.1) 
by virtue of Cardy' formula [EB]. It is in fact essential to stress that, in contrast to other 



isothermal critical amplitudes, g* is fixed by this formula, which reads [26 



c = 3nh\2 - |)2 J dVG{r) (4.2) 

where c is a central charge; G{r) is the two-spin correlation function; and ?7=0.25 is the 
anomalous scaling dimension of the spin variable a in the 2D Ising model. Furthermore, 



Eq.( [4.2|) takes on a more convenient form 

c = 377(4 - vfh'ex (4.3) 

and it should be stressed that this formula is valid at T = Tc;h ^ 0. At the Curie point the 
correlation length and the susceptibility are described by power laws 

ah) = f^h~^; x{h) = C,h^ (4.4) 

with /f , Cc being the isothermal amplitudes. On the one hand, from Eq. (p.6|) and Eq. ( |4.4|) 
it follows that 

_ 2(2-,)(4-3,) 

On the other hand, from Eq.s ( p.6|) and ( [4.3|) it is seen that the correlation length drops 
out of the product glc: 

c3: = 3.(4-,y^/,^f-5!^H-3M!) («) 

\ X J 
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Inserting Eq. (|4.4| ) into the r.h.s. of Eq. ( |4.6| ), one obtains the renormalized coupling 
constant value at the end point of the isothermal line 

g: = -i2-v)iA-3v) (4.7) 
c 

For the 2D Ising model Eq. (f4.7|) yields g* = 2737r/4. Eq. ( [4.3|) is easily seen to impose a 
constraint on the amplitudes /f and Cc 

^From Eq. ([4.8|) it follows that what we actually found, by virtue of Cardy's formula, is 
only the product Cc(/f)^. To compute these quantities separately one needs more powerful 
techniques. 



In some seminal papers p7|-p9| it was shown how to compute the isothermal amplitudes 
by making use of the Thermodynamic Bethe Ansatz and within the framework of the form- 
factor approach. In particular, in his paper Fateev obtained the following remarkable result 

m 

_ r(2/3)r(8/i5) r(i/4)r^(3/i6) _ 2270194675 u 9) 

- 4sin(7r/5)r(l/5)U7r2r(3/4)r2(13/16)J " 0-2270194675 (4.9) 



^From Eq.s (|]8D and (gj) it follows that 

Cc = 0.0731998414 (4.10) 

All these results allow one to compute exactly the two following universal combinations 
(see also ||23| and pO| ) 



Q2 = ^(4)'-^ (4.11) 



We recall that the definitions of the critical amplitudes entering Eq. ( [4.11 ) are as follows 
ll,|2|g) 



M, = B{-rf 
x(T^r, + 0) = C+r-^ 

^ = f^r-^ (4.12) 

Ms being the spontaneous magnetization. The amplitudes /f , Cc have been defined in 
Eq.([4.4|). The exact values are as follows: 6 = 15,1] = 0.25, = l,/3 = 0.125,7 = 1.75, 
B = 1.222410 ( 0), /+ = 0.567296 ( lH, see also the Appendix), and C+ = 0.962582 ( 
pT|). Substituting all these values into Eq. ( |4.11| ) yields 

Qi = 0.912648 

Q2 = 2.647714 (4.13) 
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The exact values found provide a good opportunity to test the numerical results ob- 
tained from the HT series expansions and from Monte Carlo simulations. The fair esti- 
mates obtained from the analysis of HT series in the 2D Ising model on the square lattice 
yield = 0.233 and = 0.0706 ( 0), whilst the exact results are given by Eq.s 
and ( [4.10| ). As for the universal combinations Qi^2, the series expansion analysis yields 
Qi = 0.88023, Q2 = 2.88 

Notice in conclusion that Eq.s ( [4.7] ) and ( |4.8| ) hold good also for the general case of the 
2D (7$^ 0(A^)-symmetric model for —2 < iV < 2, in particular for the minimal models of con- 
formal field theory corresponding to the discrete values of A^: N = 2 cos(— ); m = 3,A,...,oo. 
More interesting still is that certain combinations of isothermal critical amplitudes, in partic- 
ular Cc{fiY, have the same numerical values both in the pure and in the quenched disordered 
2D Ising model where c = 0.5; 77 = 0.25 [H]. 



V. CONCLUDING REMARKS 

We have proved the existence of the duality symmetry transformation d{g) in the 2D 
g^^ theory such that P{d{g)) = d\g)[3{g). Actually, this symmetry property was shown to 
result from the KW duality of the 2D lattice Ising model. 

It would be tempting but wrong to regard d{g) as a function connecting the weak-coupling 
and strong coupling regimes. For instance, this phenomenon takes place in the sinh-Gordon 
theory with f3{d{g)) = 0, the exact S-matrix being invariant under the strong- weak coupling 
duality g Sn/g ||28|. 

As a matter of fact, our proof is based on the properties of g{s),s{g) defined only for 
0<s<oo;(7*<(7<oo and therefore doesn't cover the weak-coupling region, < g < g*- 
Whether the beta-function f]{g) does have the dual symmetry in the weak-coupling region, 
remains an open question. 

In contrast to widely held views, the duality symmetry imposes only mild restrictions 
on (3{g). It means that this symmetry property fixes only even derivatives of the beta- 
function [3^'^^\g*){k = 0, 1, ...) at the fixed point, leaving the odd derivatives free. In fact, 
the duality equation d{g) = g provides yet another method for determining the fixed point, 
independently of the approach based on the equation (3{g) = 0. Another open problem is 
also that of finding a systematic approach for calculating d{g). 

In this paper the strong coupling expansion for the beta-function of the 2D Ising model 
has been developed. This approach was shown to provide reliable results for the numerical 
value of g^ which are found to be in good agreement with the best results from the HT 
series expansion. It is worth noting that g^ remains the same for the 2D Ising model with 
random bonds 0. 

The method presented in this paper may be easily extended to the 2D 0(A^)-symmetric 
g^^ theories for arbitrary A^. However, the above approach appears to be unable to provide 
reliable numerical results for the critical exponent u, at least without the combined use of 
a resummation technique. In fact, the results obtained disagree both with estimates found 
from the four-loop RG calculations, combined with the Pade-Borel resummation method 
23, and with the exact value u = A/3 given by conformal field theory P3. As a matter of 



fact, the situation with the critical exponent u looks somewhat unresolved. On the one hand. 
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all the corrections to the scaling laws in the 2D Ising model are analytical. For instance, 
the susceptibility near is given by X = C+r-^/^ + ^ On the other hand, 



corrections to scaling are known to be powers of t^'^ [^. All this would lead to a; = 1, 
in obvious contradiction to conformal field theory. Moreover, the spectrum of conformal 
dimensions of the 2D Ising model consists of just three numbers, these being (0, |, 1). That 
is the reason why the appearence of an operator with a fractional scaling dimension 4/3 
is so far unclear [^. Thus, although one can make a meaningful comparison between the 
numerical values obtained, nevertheless the apparent spread of these numerical results for 
uj does not allow one to construct a smooth interpolation for I3{g) between the two regions: 
^ < ^; and ^ > 

Regarding critical amplitudes, namely, fi,Cc,Qi,Q2, we can conclude that the known 
numerical results are in a surprisingly good agreement with the exact results. For further 
discussion, see Ref. [02l . 

We have already seen that the exact beta-function l3ising{T) is essentially non- 
perturbative. Non-perturbative terms give rise to nontrivial fixed points in contrast to the 
standard beta-functions of non-linear sigma models with continuous symmetry, in particu- 



lar for the 0(iV)-symmetric theory with N > 2 |Q. The very existence of non-perturbative 



terms is indicative of the mathematical illegitimacy of the naive analytical continuation 
from the > 2 to the N < 2 region often employed in statistical mechanics and condensed 
matter physics. 
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VII. APPENDIX 

Here we begin with giving a simple derivation of the correlation length ^ for the 2D Ising 
model from the high-temperature expansions and will then present some useful relations. 



The HT expansions up to K^(K = ^) are as follows |jT5 



X = l + 4:K + 12K^ + 104KV3 + Q2K^ + 3608irVl5 
+ 30567^75 + 4845287^7315 + 400012A:7105 (A.l) 
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^2 = 4:K + 32K^ + A88K^/3 + 20A8K^/3 + 381687^715 

+ 3946247^745 + 8994736/^7315 + 28064768/^7315 (A.2) 

Xlf^ = -2 - 32ir - 264:K^ - 4864X^/3 - 82327^^ - 553024/^715 

- 2259616/s:7l5 - 180969728/^7315 - 217858792/^7105 (A.3) 

with Xhh being the second derivative of the homogenous susceptibihty with respect to a 
magnetic field h; x ci-nd /i2 were defined in Sect. II. 

The standard RG equation for the effective temperature T is given by 

= (3lsrn,{T) (A.4) 

where ^ is the correlation length (Eq. ( |2.3| )) 



The key observation for computing ^ is to make use of the new variable s = 
exp(2/^) tanh(/C). One has to substitute Eq.s (|A.1| ) and ( |A.2| ) into Eq. (|A.5|) and then 



to rewrite the expression obtained in terms of s. At this order of approximation the proce- 
dure gives the result 



^2 = 5 + 25^ + 3s^ + 4s^ + 5s^ + 6s^ + 7s^ + 8s^ + O(s^) 



s 



1-s) 



(A.6) 



We believe the closed-form above to be an exact result (it satisfies in particular the s —>■ 
duality symmetry). However, it should not be confused with other exact formula available 
in the literature for ^, based for instance on the transfer-matrix eigenvalue approach |3| 
rather than our field-theoretic definition. 

One sees that the correlation length as given by Eq. ([A.6|) in its closed form does exhibit 
the correct behavior near the Curie point 

X T — T 2 

C = f^-; r = nn ^ = 7= = 2.269185 

^ r' T, ' ln(v^+l) 

/i+ = 1= = 0.567296 (A.7) 



where fi is the correlation length critical amplitude 

Combining Eq.s (|A.4|) and ( |A.6|) and carrying out some calculations, we arrive at a nice 



expression for the beta-function written in terms of s 



ds n / \ 2s(l - s) , . 



Substituting now Eq.s ( |A.1| ), ( |A.2| ) and ( |A.3|) into the r.h.s. of Eq. ( p.9| ) and making use of 



the variable s, one is led to some helpful relations linking s and g 
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- = s/2 - Ss^ + 39sV2 - 138s^ + 2061sV2 
9 

- 7909s^ + 122907s72 - 480012sV945 + 0(s~^) (A.9) 

s = 2/g + 24// + 264// + 2976// + 35136// + 423680// 
+ 5149824/^^ + 63275520// + 0{g~^) (A.IO) 

Rewritten in terms of T, Eq. ( |A.8| ) reads 
dT _ _ ^, 1 - 2exp(-2/T) - 2exp(-4/T) + 2exp(-6/T) + exp(-8/T) , 

raising \ ) 1 I o ^■,^-»^/' o i o ^■,^-»^/' ci ^■^,-^( Q ^ * 



^1^^ </v / 1 + 2exp(-2/T) +2exp(-6/T) -exp(-8/T) 



The exact beta-function given by Eq.s ( |A.8| ) and ( [A. Ill ) exhibits the following properties: 



(i) it vanishes just at two fixed points (FP): at the Gaussian FP s = and at the Ising 
FP s = 1; 

(ii) the derivative of the beta-function at the Ising FP gives the critical exponent of the 
correlation length ^, more precisely 

= _^^W£) _ ^ ^A.12) 
ds \s=i 

(iii) the beta-function is covariant under the duality transformation, namely 

Pisingis) = -S^Pi.ingi^) (A. 13) 

and this equation has been referred to as the consistency condition for the beta-function 

(iv) these properties fully agree with results from the HT series expansions. 

The duality property (iii) implies that the RG equation under discussion takes the same 
form both in the dual and in the original variables. We note in passing that the well-known 
self-dual beta-function obtained within the Migdal approximation for d = 2, which reads 



1 35 



2 1 

pMigdai (T) = -T - — sinh — In tanh — (A. 14) 



does not satisfy either properties (ii) and (iv) p5| . 

To end with, we prove some useful relations concerning the duality transformation g = 
d{g) introduced in Sect. II. 

i) Let us show that d'{g*) = ±1. First of all, from the definition 

dig)^s-\^^) (A.15) 
we have that d{d{g)) = g. Differentiating this with respect to g one obtains 
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d{d{g)) = g- dig*) = g* 

d\d{g))di{g) = 1; g = g* 

df{g*)df{g*) = 1 (A.16) 



showing that indeed df{g*) = ±1. 




(A. 17) 



and at the fixed point we arrive at {d{g*) = g*): 



d"{g*){d{g*)f + d\g*)d"{g*)=Q 



(A.18) 



If d'{g*) = —1 we have an identity d"{g*) = d"{g*). In the opposite case from d'{g*) = +1 
it follows that d"{g*) = 0. Proceeding in the same way, it is easy to see that all higher 
derivatives vanish identicatically at the fixed point. This, alone, does not imply that d{g) = 
g, since it could be that d{g) = g + f{g) where f{g) is some nonanalytic function having 
vanishing derivatives at g*. Let us then assume that f{g) is the first term of an asymptotic 
expansion of d{g) around the fixed point, so that f{g) ^ as g g*. Then remembering 
that d{d{g)) = g, 



g = d{d{g)) = d{g + f{g)) = g + f{g) + f{g + f{g)) ^ g + 2f{g) 



(A.19) 



and we do obtain f{g) = 0. 
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